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Abstract 



In supersymmetric models, the CP asymmetry produced in the decay of the 

lightest right-handed neutrino, = e, can be written as a function of weak scale 

0^ . parameters. We introduce a way of separating e into contributions from the 

various weak-scale phases, and study the contribution of potentially measurable 

neutrino phases to leptogenesis. We find that the Majorana phase </>', which 

could have observable effects on neutrinoless double beta decay, is important for 

e unless there are cancellations among phases. If the phase 5 can be measured at 

\Q ' a neutrino factory, then it contributes significantly to e over much of parameter 

O . space. 
(N 
O 

23 
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43 ! 1 Introduction 



After the discovery of neutrino oscillations[l,2], leptogenesis [3] stands as one of the 
most appealing explanations for the observed Baryon Asymmetry of the Universe 
(BAU) [4]. One of the crucial ingredients [5] for this mechanism is CP violation in 
the leptonic sector. However, there is no indication for it so far, although it could 
perhaps be seen at a neutrino factory or in neutrinoless double beta decay. It is there- 
fore interesting to investigate whether there is any relation between the CP violation 
required for leptogenesis and the phases that could be measured at low energies in the 
neutrino sector. 

The leptogenesis scenario relies on the seesaw model [6] for neutrino masses, that 
is usually analyzed in terms of high-energy parameters, not accessible to experiments. 
So, the resulting predictions are (texture) model-dependent. The above question has 
been addressed in such an approach [7-12]. Instead, we parametrize the seesaw in 
terms of weak scale variables [13]. This gives us a model-independent formulation of 
leptogenesis in terms of low energy inputs, in which we can study the above question. 

The aim of this paper is to quantify, in a model independent way, the relation of 
the CP violation required for leptogenesis to the measurable low-energy phases. We 
express the CP asymmetry of leptogenesis as a function of real parameters and phases 
at the weak scale, and then introduce a definition of "phase overlap" between the 
leptogenesis phase and the individual low energy phases. This definition is not the 
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only possible one, but has linearity properties and is calculable. It is motivated by the 
notion of vector space, spanned by low energy phases ("basis vectors"), in which the 
CP asymmetry of leptogenesis is a "vector" . The relative importance of a low energy 
phase for leptogenesis would then be the "inner product" of the leptogenesis "vector" 
with the relevant "basis vector" . We will not be able to construct such a vector space, 
but it is a useful analogy to keep in mind. 

The paper is organized as follows: the next section introduces CP violation in the 
leptonic sector. Section 3 contains the basic concepts of the supersymmetric see-saw 
and the generation of the BAU by leptogenesis, from a top-down point of view. In Sec- 
tion 4 we review the procedure to reformulate the see-saw mechanism from a bottom-up 
perspective. This will allow us to study leptogenesis in terms of low energy data, open- 
ing the possibility of relating, in a straight-forward way, the Baryon Asymmetry of 
the Universe with the CP violation measurable at neutrino factories. In section 5 we 
develop the general formalism to study quantitatively the above-mentioned relation. In 
Section 6 and 7 we show the results of our analysis, first for a particular case, and then 
for a more general case. In section 8 we present a self-contained summary and conclu- 
sions. Finally, we include an appendix with the procedure to evaluate numerically the 
contributions from the low-energy phases to leptogenesis. 

2 Flavour and CP violation in the leptonic sector 

In the last few years, the Superkamiokande collaboration [2] has provided compelling 
evidence that neutrinos have mass and oscillate. More recently, the SNO collabora- 
tion [14] has confirmed the oscillation hypothesis, and the first neutral current data 
[15] seem to favour the large angle MSW (LAMSW) solution to the solar neutrino 
problem [16]. These results, combined with those from a series of other experiments 
[17], have allowed to measure fairly well the mass splittings and mixing angles relevant 
for solar and atmospheric neutrino oscillations. In addition to this, other experiments 
have provided bounds on neutrino parameters from electron antineutrino disappearance 
(CHOOZ)[18], the non-observation of neutrinoless double beta decay [19], the shape 
of the tritium beta decay spectrum [20], and different cosmological and astrophysical 
considerations. However, no evidence has been found so far for CP violation in the 
leptonic sector. 

The search for leptonic CP violation is theoretically motivated by several facts. 
First, the discovery of CP violation in the leptonic sector could shed some light on 
the mechanism that generates neutrino masses and perhaps hint at some underlying 
structure. Secondly, the observation of CP violation in the quark sector, (in the neutral 
kaon system, e'/e, and in B — > ipK s ), encourages the search for CP violation in the 
neutrino sector. If there exists a symmetry relating quarks and leptons, these exper- 
imental results would point to CP violation also in the leptonic sector. Furthermore, 
particular models would give definite predictions that could be contrasted in the fu- 
ture. Lastly, and most importantly for the purposes of this paper, CP violation in the 
leptonic sector could be related to the observed Baryon Asymmetry of the Universe. 
This is possible in the context of the see-saw mechanism. 
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On the experimental side, the leptonic version of the CKM phase can be detected 
by comparing transition probabilities for neutrinos and antineutrinos: 

A _ P{y a -> vj) - P{y a -> vj) (1) 

P{v a -> v p ) + P(u a -> Up) 

Unfortunately, it is not possible to measure such an asymmetry with the natural sources 
of neutrinos, i.e. the Sun and pions decaying in the atmosphere, since the "beam" 
cannot be switched from v to v. Hence, a lot of effort is being bestowed on the 
design of a neutrino factory [21,22]: an intense muon source to produce a high-intensity 
neutrino beam. In the muon storage ring, muons decay to produce muon neutrinos 
and electron antineutrinos. Whereas a muon neutrino would produce a muon in the 
detector, the oscillation of an electron antineutrino to a muon antineutrino would 
produce an antimuon. This antimuon (a "wrong sign" muon) would be a clear signature 
for oscillation, and P(P e — > v^) could be determined. One of the advantages of a 
neutrino factory is that the muons in the storage ring can be replaced by antimuons. 
This makes possible the measurement of P[y e — > and hence the CP asymmetry. 
In practice, detecting CP violation in the neutrino sector is not an easy task [23,22], 
since the beam has to go through the Earth, that is CP asymmetric. In consequence, 
the matter effects on the oscillation pattern can obscure the CP violation intrinsic to 
neutrinos. 

If neutrinos have Majorana masses, as predicted by the seesaw mechanism, there 
are also "Majorana" phases, in addition to the "Dirac" phase that could be detected at 
a neutrino factory. Neutrinoless double beta decay could be sensitive to these phases 
(see however [24]). This lepton number violating, but CP conserving, process probes 
the Majorana neutrino mass matrix element between v e and is e , which depends on the 
masses and mixing angles, and also on the Majorana phases. Neutrinoless double beta 
decay is not observed at the moment. However, experiments which should see a signal, 
for the currently favored masses and mixing angles (LMA), are being discussed [25]. 

The see-saw mechanism [6] consists on adding three right-handed neutrinos to the 
Standard Model (SM) particle content, singlets with respect to the SM gauge group, 
and coupled to the Higgs doublet through a Yukawa coupling. Then, a Majorana mass 
term for the right-handed neutrinos is not forbidden by the gauge symmetry, and can be 
naturally much larger than the scale of electroweak symmetry breaking. These simple 
assumptions are enough to produce neutrino masses naturally small 1 . Furthermore, if 
CP is violated in the leptonic sector, the decay of the right-handed neutrinos in the early 
Universe produces a lepton asymmetry [3,28] that will be eventually reprocessed into 
a baryon asymmetry by sphalerons [29] . This leptogenesis scenario will be discussed in 
more detail in section 3. 

In supersymmetric models, the seesaw mechanism can induce flavour violating pro- 
cesses involving charged leptons that could be observed in the future [26]. The neutrino 
Yukawa couplings generate off-diagonal elements in the slepton mass matrix, via renor- 
malization group running. These flavour violating mass terms contribute inside loops 

Nevertheless, this minimal model has a serious hierarchy problem: the right-handed neutrinos 
produce a (large) quadratically divergent radiative correction to the Higgs mass. Therefore, in what 
follows, we will restrict ourselves to the supersymmetric version of the see-saw mechanism. 
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to processes such as /x — > e7, r — > 7/7 and r — > e7. This has been extensively studied 
from various theoretical [30,12] and phenomenological [31,32,27] perspectives. The cur- 
rent experimental bound [33] on \x — > ey imposes some restrictions on the parameter 
space of the SUSY seesaw. It is anticipated that the sensitivity to r — > /ry and fi — > &-/ 
could improve by as much as three orders of magnitude [34] in forthcoming years. This 
would provide interesting information about the flavour structure of the SUSY seesaw, 
irrespective of whether lepton flavour violation is observed or not. 

In this paper, we will concentrate on the possible relation of the CP asymmetry in 
the leptonic sector with the Baryon Asymmetry of the Universe, in the framework of 
the super symmetric leptogenesis. We suppose that the BAU is generated in the out- 
of-equilibrium decay of the lightest right-handed neutrino. The CP violation that gives 
rise to the BAU is not straight-forwardly related to the CP violation that could be 
observed at low energy. This has been carefully and elegantly discussed, using Jarlskog 
invariants in [7] . These authors have also studied the relation of the leptogenesis phase 
to low energy phases for specific high scale models and various solar solutions [8] . They 
have an analytic approximation similar to ours, and our results seem to agree where 
they overlap. Correlations between leptogenesis and low energy parameters have been 
studied in left-right symmetric models where the Yukawa couplings are small [10], and 
in various Grand Unified Theories [11]. In certain classes of top-down models, it has 
been found [12] that the leptogenesis phase is unrelated to the MNS phases; these 
textures correspond to the third case we study, in section 7.2, so the conclusion [12] 
that 5 makes little contribution to leptogenesis agrees with our result. The goal of this 
paper is to investigate the interplay between the CP violation at very high energies and 
at low energies, in a model independent way. We will also comment on the prospects 
to observe CP violation at a neutrino factory or in neutrinoless double beta decay, in 
view of the measured BAU, and inversely, what could be inferred about the BAU if 
CP violation is observed at low energy. 

3 The see-saw mechanism and leptogenesis: the 
top-down approach 

The supersymmetric version of the see-saw mechanism has a leptonic superpotential 
that reads 

W lep = e c R T Y e L ■ H d + v c R T Y v L ■ H u - \v R T Mv R , (2) 

where Lj and em (i = e, //, r) are the left-handed lepton doublet and the right-handed 
charged-lepton singlet, respectively, and Ha (H u ) is the hypercharge —1/2 (+1/2) Higgs 
doublet. Y e and Y u are the Yukawa couplings that give masses to the charged leptons 
and generate the neutrino Dirac mass, and M. is a 3 x 3 Majorana mass matrix that 
does not break the SM gauge symmetry. We do not make any assumptions about the 
structure of the matrices in eq.(2), but consider the most general case. Then, it can be 
proved that the number of independent physical parameters is 21: 15 real parameters 
and 6 phases [36]. 
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It is natural to assume that the overall scale of M., denoted by M, is much larger 
than the electroweak scale or any soft mass. Therefore, at low energies the right-handed 
neutrinos are decoupled and the corresponding effective Lagrangian reads 

5C lep = e c R T Y e L ■ H d - \^Y V L ■ H U ) T M.-\Y V L ■ H u ) + h.c. (3) 

So, after the electroweak symmetry breaking, the left-handed neutrinos acquire a Ma- 
jor ana mass, given by 

M v = m D T A^ _1 m D , (4) 

suppressed with respect to the typical fermion masses by the inverse power of the large 
scale M. 

We will find convenient to work in the flavour basis where the charged- lept on 
Yukawa matrix, Y e , and the gauge interactions are flavour-diagonal. In this basis, 
the neutrino mass matrix, A4 U , can be diagonalized by the MNS [37] matrix U, defined 
by 

U T M V U = di&g(m 1/1 ,m 1/2 ,m V3 ) = D Mv , (5) 
where U is a unitary matrix that relates flavour to mass eigenstates 



(6) 



and the m Ui can be chosen real and positive. Also, we label the masses in such a way 
that m Ul < m U2 < m U3 . We will assume throughout the paper that the light neutrinos 
have a hierarchical spectrum. We do not consider the inverse hierarchy, which may be 
more difficult to match with the neutrinos detected from SN1987A [39]; we anticipate 
that the contribution of 5 to the leptogenesis phase could be suppressed by powers 
of m Ul in this case. The CP asymmetry required for leptogenesis is suppressed for 
degenerate light neutrinos [38], so we neglect this possibility. 
U can be written as 

U = V-diag(e-^ 2 ,e-^' /2 ,l) , (7) 

where and 0' are CP violating phases (if different from or n) and V has the ordinary 
form of the CKM matrix 















\Vt) 







/ C13C12 C13S12 s 13 e~ lS \ 

V = -C23S12 - S23Sl3Cl 2 e l<5 C23C12 - S 2 3Sl3Sl2e l5 S23C13 

V S23S12 - c 2 3S 13 c 12 e tS -S23C12 - c 23 si3Si2e l<5 c 23 ci 3 / 



(8) 



It is interesting to note that the neutrino mass matrix, Ai u , depends on 9 parameters: 
6 real parameters and 3 phases. Comparing with the complete theory, we discover 
that some information has been "lost" in the decoupling process, to be precise, 6 real 
parameters and three phases. We will return to this important issue later on. 

Another remarkable feature of the see-saw mechanism is that it provides a natural 
framework to generate the baryon asymmetry of the Universe, defined by i] B = (ub — 



5 



n§) I s, where s is the entropy density. This quantity is strongly constrained by Big Bang 
Nucleosynthesis to lie in the range rj B ^ (0.3 — 0.9) x 10~ 10 , to successfully reproduce 
the observed abundances of the light nuclei D, 3 He, 4 He and 7 Li [40]. As was shown by 
Sakharov, generating a baryon asymmetry requires baryon number violation, C and CP 
violation, and a deviation from thermal equilibrium. These three conditions are fulfilled 
in the out-of-equilibrium decay of the right-handed neutrinos and sneutrinos in the 
early Universe. For conciseness, and since we are concerned only with supersymmetric 
leptogenesis, in what follows we will use right-handed neutrinos, and the shorthand 
notation v R , to refer both to right-handed neutrinos and right-handed sneutrinos. 

Let us briefly review the mechanism of generation of the BAU through leptogenesis 
[3,28]. At the end of inflation, a certain number density of right-handed neutrinos, n UR) 
is produced, that depends on the cosmological scenario. These right-handed neutrinos 
decay, with a decay rate that reads, at tree level, 

r A = r(u Ri - UE U ) + Y{v Ri - UK) = ^(Y u Y u %M t . (9) 

The out of equilibrium decay of a right-handed neutrino v Ri creates a lepton asymmetry 
given by 

n t ~ n i n UR + n„ R 
Vl = = €i K. (10) 

s s 

The value of (n UR + no R )/s depends on the particular mechanism to generate the right- 
handed (s) neutrinos. On the other hand, the CP- violating parameter 

f , = ^-|^, (id 

where F D . is the CP conjugated version of T D ., is determined by the particle physics 
model that gives the masses and couplings of the v R . Finally, k is the fraction of the 
produced asymmetry that survives washout by lepton number violating interactions 
after v R decay. To ensure k ~ 1, lepton number violating interactions (decays, inverse 
decays and scatterings) must be out of equilibrium when the right-handed neutrinos 
decay. In the case of the lightest right-handed neutrino u Rl , this corresponds approx- 
imately to T Dl < H\ T ~ Ml , where H is the Hubble parameter at the temperature T, 
and can be expressed in terms of an effective light neutrino mass [28,41], mi, as 



^g^r Dl = {Y V Y J) J^jL < 5 x 10- 3 eV. (12) 



This requirement has been carefully studied [28,41,50]; the precise numerical bound on 
mi depends on Mi, and can be found in [41]. 

The last step is the transformation of the lepton asymmetry into a baryon asym- 
metry by non-perturbative B+L violating (sphaleron) processes [29], giving 

C 

Vb = ^T—[Vl, (13) 

where C is a number 0(1), that in the Minimal Supersymmetric Standard Model takes 
the value C = 8/23. 
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In this paper, we assume that a sufficient number of vr were produced — thermally, 
or in the decay of the inflaton, or as a scalar condensate of vrs, or by some other 
mechanism. We will concentrate on the step of leptogenesis that is most directly 
related to neutrino physics, namely the generation of a CP asymmetry, e, in the decay 
of the right handed neutrinos. It is convenient to work in a basis of right-handed 
neutrinos where M. is diagonal 

M = diag(M 1 ,M 2 ,M 3 ), (14) 

with Mi real and < Mi < M 2 < M 3 . In this basis, the CP asymmetry produced in 
the decay of u Ri reads 



C 7. 



8tt [YyYjljj 



ElHi{[Y,Y,t]j}/^, (15) 



where [43] 



/(a;) = v /i(^- + ln[i±^). (16) 
\x — 1 L x J / 



1 +x 



Here, we will assume that the masses of the right-handed neutrinos are hierarchical. We 
also assume that the lepton asymmetry is generated in the decay of the lightest right- 
handed neutrino. This second assumption is critical; if the asymmetry was generated 
by the decay of vr 2 or z/r 3 , it would depend on a different combination of phases. This 
assumption is also dubious if the v R are produced thermally, because the u Rl mass, 
Mi, is severely constrained in SUSY models. To get a large enough baryon asymmetry, 
Mi > 10 s GeV is required [38,42], but Mi must be less than or of order the reheat 
temperature T re h. To avoid overproducing gravitons in the early Universe, T re h is 
required to be < 10 9 - 10 10 GeV [44]. 

With these approximations, the CP asymmetry is 

^-lw^V lm{lY ^' } i^ (17) 



3 Mi 



j—lm{[Y u MjY/] u }. (18) 



8n(H^ [Y„Y U % 

The CP asymmetry depends on quantities that appear in the superpotential of 
the complete theory, eq.(2), and that are not directly measurable with experiments. 
However, these quantities can be related to neutrino and sneutrino parameters, as we 
will discuss in the next section. One of the goals of this paper is to implement in an 
explicit way these constraints on the CP asymmetry, eq.(18). 



4 The see-saw mechanism and leptogenesis: the 
bottom-up approach 

Our starting point will be the procedure presented in [13]. In the basis defined in 
section 3, where the charged lepton mass matrix and the right-handed Majorana mass 
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matrix are diagonal, the neutrino Yukawa coupling must be necessarily non-diagonal. 
However, it can be diagonalized by two unitary transformations: 

Y„ = V^D Y V L . (19) 

It is clear that the CP asymmetry depends just on Vr and Dm- These quantities 
are related to the physics of the right-handed neutrinos and are not directly testable 
by experiments, since they are related to very high energy physics. However, there 
is a reminiscence of Vr and D M in the low energy neutrino mass matrix that can be 
exploited to obtain information about the high-energy physics from the neutrino data. 
Substituting eq.(19) in the see-saw formula, M u = m D T A / f _1 m D , one obtains 

Dy'Vt^VlDy 1 = V2£tf Vj = M~\ (20) 

From this equation we can solve for V R and D M in terms of M u , D Y and V L . M u 
is constrained by neutrino experiments, whereas D Y = diag(yi, y 2 , yz) and Vl are 
unknown parameters at this stage. We choose a parametrization of the unitary matrix 
Vl such that 



Vr 



/ r L „L „L Q L „—i(pi2 a L „—i(pi3 

c 13 c 12 c 13 6 12 e 6 13 e 

_ r L Q L p iipi2 _ Q L q L „L p i{tprj,-tp2z) r .L „L _ „L „L „L p -i(ipi2-<Pi3+</>23) Q L „L p -i<pn 

°23*12 e i, 23 6 13 c 12 e c 23 c 12 i, 23 6 13 i '12 e i '23 c 13 e 



„L „L <3 i(ipi2+</ 3 23) n L Q L n L „iipi3 a L „L „iip23 n L Q L Q L ^{ipvA^^u) r .L r .L 

\ i, 23 i '12 e c 23' b 13 c 12 e i '23 c 12 e c 23 i '13 i '12 e c 23 c 13 

(21) 

where cfj = cos 9fj and sfj = sin^-, being 9^ the angles in the Vl matrix. 

In certain scenarios, the parameters Dy and Vl can be constrained experimentally. 
For example, in a scenario of minimal SUGRA, with just the MSSM+3z/rs below the 
GUT scale, D Y and Vl can in principle be extracted from the radiative corrections 
to the left-handed slepton mass matrix, since the corresponding RGE depends on the 
combination Y^Y^ = vId y Vl- To be more precise, at low energies the left-handed 
slepton mass matrix reads, in the leading-log approximation 

Ks) y - (diagonal part) £> + ^(3m^ + A^Y^Y^ log . (22) 

The off-diagonal terms in m|. manifest themselves in processes like \i — > or r — > /ry, 
that could be observed in the near future 2 . In addition to this, at tree level the 
three sneutrino masses are degenerate. However, radiative corrections induce a non- 
universality among the masses that could perhaps be measured experimentally [46]. All 
these measurements could be used to disentangle some information about the neutrino 
Yukawa matrix and the right handed masses from radiative corrections. See [47] for a 
recent analysis of £j — > in this approach. 

For leptogenesis we are particularly interested in the phases of Y u , that are in turn 
related to the phases in the left-handed slepton mass matrix. It is then an important 

2 The upper bounds on the off-diagonal entries of Y^Y^ log ^f- in the scenario of mSUGRA 
with the MSSM+3i/RS below the GUT scale apply for a wide class of models, since one does not 
expect cancellations among the different terms in the RGEs, or with the off-diagonal elements of the 
tree-level slepton mass matrix. 
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issue to measure the phases in m|_. The electric dipole moments (EDMs) of the 
electron and muon are CP violating observables that could provide information about 
these phases. However, this CP violation is flavour conserving and could come from 
another flavour conserving sector of the theory, like the charginos and neutralinos, 
instead of the Yukawa couplings. Furthermore, the EDMs care about the relative 
phases between the charged leptons and the sleptons, of which the see-saw only induces 
one. So, the contribution induced by the see-saw would be suppressed by small angles 
and Yukawa couplings [48]. Therefore, to constrain the phases in the Yukawa couplings 
with the EDMs, one has to make certain assumptions about the soft SUSY breaking 
lagrangian. A more detailed discussion of obtaining information about the complete 
theory from low energy data can be found in [13]. 

The V L , D Y low-energy parametrization has several advantages. If we treat the 9 
parameters of the neutrino mass matrix as "known", there are 9 remaining unknown 
variables in the seesaw: three phases and six real numbers. Possible parametrizations 
of these unknowns are Dy and Vl, Dj^ and the orthogonal complex matrix R [27], or 
as in [47] . The angles and phases of Vl are related in a simple way to the lepton flavour 
violating slepton mass matrix entries. These off-diagonal (in the charged lepton mass 
eigenstate basis) entries are currently constrained and could possibly be determined by 
radiative lepton decays £j — > ^7. The eigenvalues of Dy are more difficult to determine 
experimentally. However, we do measure the Yukawa matrix eigenvalues for the quarks 
and charged leptons, so we can make theoretical guesses of the Y v eigenvalues with more 
confidence than e.g. guessing the v R Majorana masses. 

It is convenient for our leptogenesis analysis to parametrize the sneutrino mass 
matrix with Dy and Vl- It would be more correct to express the lepton asymmetry in 
terms of the magnitude and phases of slepton mass matrix elements 3 . Alternatively, 
there is an intermediate parametrization, which can be useful for analytic estimates. 
The parameters we use, V L and D Y) determine Y\Y U rather than Yj.^Y^ log ^ , 
which is the expression that appears at leading log. It is easy, though, to relate Vl and 
Dy to YlY Ukj log 

m G ut ' Noting that 



Y ^ k bg i^ Y ^ = (Y ^ )y ' = ^ D ^ Vl ^ (23) 
(rUT M m , 



where Y ki = Y ki ^J\og^^ and M k = M fc log-^^, it is possible to rewrite eq.(20) 
but using tilded parameters. So, one could parametrize the see-saw mechanism with 
the neutrino mass matrix, Ai v , and Vl, Dy, that are directly related to the leading-log 
approximate solution of the left-handed slepton RGEs. Also, from the definitions, it 
is straightforward to relate Vl and Dy with their tilded-counterparts. However, since 
SUSY has not yet been discovered, we use Vl and D Y , with the knowledge that we can 
calculate [m~] from these parameters. This choice will be important when we discuss 
phase overlaps. 



3 We will follow this approach in a subsequent publication [49] . 
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We turn now to expressing the CP asymmetry in terms of neutrino masses, the 
MNS matrix, and other unknown parameters encoded in D Y and V L . We can make 
an analytic approximation indicating the dependence of the CP asymmetry e on our 
low energy parameters. To derive these estimates, we first assume M 3 ^> M 1 and 
V\ ^ 2/2,2/3- Then we assume that [.M^-M -1 ]!! is the largest element of M-^M.^ 1 , 
in the basis where Y v is diagonal. As we will see, this is usually reasonable. 

If a matrix A has a zero eigenvalue, then the remaining two eigenvalues are 

Ai, \ 2 = \ {TrA ± \J (TrA) 2 - 4 (A n trA + detA - |A 12 | 2 - |A 13 | 2 )} , (24) 

where Tr is the trace of the 3-d matrix, and tr and det are defined on the 2-3 subspace. 
In the limit where M3 — > 00, this formula can be applied to the hermitian matrix 

M M = ^ = - { . (25) 

To obtain simple expressions, we would like to expand eq. (24) in small dimensionless 
parameters. So to avoid confusion, we scale a factor y\ out of Mr x M.. 

The largest eigenvalue of A will be of order m 2 , as can be seen by defining 77 = 
ytDy 1 = dia^ {771, 7/2, and 

A = VlM v Vl = V[WD Mv U^Vl = W*D Mv W* , (26) 



which gives 



(#0)4 • W 



We take 7/3 = 1. The matrix W = VjJJ is the rotation from the basis where the vl 
masses are diagonal to the basis where the neutrino Yukawa matrix Yj Y„ is diagonal. 

The dominant contributions to the matrix elements of A can be calculated as an 
expansion in 772 and 773. Only An is zeroth order in 772 and 773, so generically An ^> 
A22, A 33 . Then from eq. (24), the lightest RH neutrino has mass 



(Mi) 2 ~ , (28) 



A 
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and the associated eigenvector will be 




A* 
^11 



v x ~ A21 x — = 77 2 A* 12 x — . (29) 
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We can use this eigenvector to evaluate eq. (18), and find 

„ 3yfAfi lm l [AD Y AW Y A T ] n \ = 3y| lm \ ^ W ^ m 

8 7 r[AL>2 A ] 11 \ [A77At77A r ] n J 8tt 1 W^m* \ En W? n m 



Vn 



(30) 
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where we have dropped terms of order 772 and 773, and recall that W is the rotation 
from the basis where the v L masses are diagonal to the basis where YjY„ is diagonal. 

It is important to notice that the CP asymmetry depends only on the first row of 
the matrix W, that in turn depends only on the first row of Vl. In the parametriza- 
tion that we have chosen for Vl, eq.(21), the first row depends on two angles and 
two phases. Therefore, at the end of the day, the CP asymmetry depends on the 
neutrino mass matrix and five unknown parameters: y±, two angles and two phases. 
Note that for generic A, the order of magnitude of e is fixed by y\. For the GUT- 
inspired value yi ~ m u /m t ~ 10 -4 , e ~ 1CT 9 unless there is some amplification in 
Im{[AA t A] n A* 1 }/ QAA^nlAnl 2 ). 

5 Phases for leptogenesis 

From the previous discussion, we find that in the parametrization we have chosen, the 
CP asymmetry depends on five phases, namely the phases in the MNS matrix, 5, <fi 
and </>', and the phases in the first row of the Vl matrix, </?i2 and </?i3. In this section we 
would like to study the relative importance of these phases on the CP asymmetry, and 
whether any of them could be considered as the "leptogenesis phase", i.e. the phase 
that is fully responsible of the CP asymmetry. 

To this end, we first introduce a definition of "overlap" between the "leptogenesis 
phase" and the low energy phases. At the end of the section, we will discuss issues 
raised by our definition. 

We define the contribution to the CP asymmetry from a phase a (this is not quite 
what we call a phase overlap) such that the total CP asymmetry is the sum of the 
different contributions: 

e = e<5 + e</> + £</>' + + e ^i 3 • (31) 

To obtain a decomposition of the CP asymmetry in this way, and give a more 
precise definition of the different contributions, we Fourier expand the CP asymmetry: 

e = Ajkimn sin(j<5 + k<p + 1(f)' + rrupu + mp 13 ) . (32) 

j,k,l,m,n 

This summation can be split in 

(33) 

a a</3 a</3<7 a</3<~f<p a</3<~f<p<a 

with {a, (3, 7, p, a} elements of the ordered set {5, <p, </>', </?i2, ^13}- Note that the subindices 
of the C's are ordered, so 5 < <fi < 0' < </?i2 < y?i3. (to avoid double counting, only 
exists, and C^s does not.) Some of the terms in the summation are: 

C 5 = J2 4;oooo sm(jS) (34) 
C54, = Y A Jkooo sin(j5 + k4>) (35) 
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C s <w = 12 A jkioo sin(j(5 + k<p + l<f>') 



(36) 



and so on. 

We can now rewrite the summation eq.(33) in a way that resembles eq.(31). It is 
clear that C$ is a contribution from S to the CP asymmetry, so it must be one of the 
terms in 6$. On the other hand, is a contribution from S, but also from 0, and it is 
not possible to conclude whether it is a contribution from 5 or from <fi. So, we will say 
that Cstj, contributes in to and in to This rationale can be applied 

to the rest of the terms in the expansion eq.(33) to finally obtain 

111 1 

e«5 = C 6 + - C sp + o 12 C sp-r + 7 12 C Wip + 7 12 c sp 1P a , (37) 

P /3<7 /3<7<P /3<7<p<cr 

and similarly for €</>, €</>/, and e^. It can be checked that with this decomposition, 
eq.(31) holds. 

In this analysis, we are only concerned with the relative contributions of the dif- 
ferent phases to the CP asymmetry, and not with the overall magnitude of the CP 
asymmetry itself (that is essentially determined by the unknown parameter yi). So, we 
normalize the different contributions to 1, and define "phase overlap" as the normalized 
contribution from a phase to the CP asymmetry. This quantity measures the relative 
importance of that phase for the CP asymmetry compared to the rest of the phases. 
For instance, the "overlap" of 5 with the leptogenesis phase is: 

S = -M= , (38) 

which satisfies EO„ = 1 and < O a < 1. Had we chosen a linear normalization, i.e. 
O a = e a /e, in some regions of the parameter space \O a \ could be larger than one, due to 
cancellations among the different e a 's. So, we prefer to use a quadratic normalization, 
that satisfies < O a < 1 to better represent the fractional contribution of the phase a 
to e. 

The overlap defined in eq.(38) measures the importance of 5 for e, provided that the 
phases in the expansion are independent and "orthogonal" . In any parametrization of 
the seesaw, six phases are required, so independence is automatic. The importance of 
"orthogonality" can be understood by analogy with linear algebra, where a vector can 
be uniquely decomposed in its components along a given orthonormal basis. Similarly, 
the phases of our parametrization must be "orthogonal" , as well as independent, to have 
a unique definition of the fraction of e due to 5. However, a mathematical definition 
of "orthogonality" is difficult, perhaps impossible, because we do not have an inner 
product between phases. So we opt for a physical notion: we assume as "orthogonal" 
the so-called "physical" phases, the phases that could be measured at low energy — in 
practice, or in principle in the best of all physicists worlds. Notice that the choice of 
low energy phases is important — had we parametrized with the phases of U and W, 
then from eq. (30), e depends only on the phases of Wu and is independent of the MNS 
phases. Similarly, if the seesaw is parametrized using U and the complex orthogonal 
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matrix R, The MNS matrix cancels out of the equation for e, and the 5 dependence of 
e is buried in R. 

The measurable phases of the slepton sector are those of the slepton mass matrices, 
so we should expand e on 5, 0', 0, and the phases of [m~J. However, we find it convenient 
to parametrize the seesaw in terms of Vl and D Y , rather than [m|], as discussed in 
section 4. The [m|] are therefore functions of the real angles of Vl, as well as the 
phases of Vl, so ipu and </?i3 are not quite the correct physical phases. We expect this 
choice to have little effect on the relative importance of 5 and 0' for leptogenesis. Vl 
is closely related to the slepton mass matrix, and in the limit that the real angles in 
Vl are small, the phases of [m|] are those of Vl (or Vl) in leading log. 

Finally, this definition of overlap is statistical . 0$ gives the probable importance 
of 5 for e, assuming all phases are unknown and 0(1). But if, for instance, 5 — 0, Og 
will be non-zero, because parts of the crossed terms {C$ a , Cs^p, •••} are included in 0$, 
and these terms could be non-zero due to the other phases. 

To avoid possible confusion, observe that our notion of "leptogenesis phase" differs 
from the one introduced in [42], who write e = e max 5 e ff where e max is the upper bound 
on e. The asymmetry is the imaginary part of a complex number e = Im{|e c |e* p } = 
I e c I ship, so we interpret the "leptogenesis phase" to be p. 

The definition, eq. (38), of the fraction of e that is due to 5, depends on eight 
unknowns: three real angles and the five phases. We assume that the real angles could 
be measured, so we present results for different fixed values of the angles. We take 
random values of the phases, linearly distributed between and 2n, and make scatter 
plots of the overlaps {O a }. If most of the points are distributed at \O a \ 2 > .3, we 
conclude that the phase a contributes significantly to e (a here represents any phase 
among {5, 0', 0, ipi 2 , fu})- Notice however that this is a statistical statement, based 
on choosing all the phases randomly and large. 

6 The case V L = 1 

In this section we particularize the previous study to the case Vl = 1. In this case, 
there is no flavour or CP violation induced radiatively by right-handed neutrinos in 
the slepton mass matrices. Since we have fixed the Vl matrix, the number of unknown 
parameters is reduced, and the CP asymmetry depends just on the neutrino mass 
matrix and y\, the lightest eigenvalue of Y^Y^. 

When Vl = 1, only the phases in the MNS matrix (5, <fi and 0') are relevant, hence 
the analysis of the previous section is greatly simplified. The CP asymmetry can be 
written as the sum of the contributions from the phases 5, and 0', 

e = e s + e<p + . (39) 

As in the previous section, we Fourier expand e, yielding 

e = A 3ki sin (i<5 + H + 1(f)') = C 5 + C + CV + C 54> + CV + C H , + C w , (40) 

3,k,l 
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with C a , C aj 3 and Cs^ as in eqs.(34)-(36). Then, the contribution from the phases 5, 
and <$' to the CP asymmetry read 

e<5 = Cs + -^{Cscf, + Cs<j>>) + -Cs^' 

1 1 

£<P — C(/> + -^{Cs(t> + C^i) + —Csc/xt)' (41) 

e^' = C^' + -(Csc/)' + C^/) + -Cs^' ■ 

As before, and since we are only concerned with the relative contributions from the 
different phases to e and not with the overall magnitude, we define the "phase overlaps" 
as: 

O a = . (42) 

V e l + 4 + e l> 

where a = 5, 0, 0'. With these definitions, the following identity holds: 

{Osf + m 2 + (£V) 2 = 1 • (43) 

In Figure 1 we show the numerical results for different CHOOZ angles. We show the 
results for the LAMSW solution to the solar neutrino problem and the mass hierarchy 
m Ul : m U2 : m„ 3 = 1CT 2 : 0.1 : 1. Each point corresponds to a random value of 
the phases 5, and 0' between and 2n. In view of eq.(43), we find convenient to 
present the results on a triangular plot, where the distance to the sides of the triangle 
corresponds to the "phase overlaps" squared, defined in eq.(42) (see upper left plot). 

When the CHOOZ angle is close to the experimental bound (upper right plot) over 
most of the parameter space the relevant phases are 5 and 0', and their contributions 
are approximately equal. In this case, the phase is essentially irrelevant, except for 
a few points that correspond to 25 — 0' ~ 0, n. On the other hand, when the CHOOZ 
angle is moderately small (lower left plot), we find points scattered over the whole 
triangle: the three phases are relevant in this case. One can also see from the figure 
that the points seem to follow a circular pattern. We will come back to this issue 
later on. Finally, when the CHOOZ angle is very small (lower right plot), the relevant 
phases are and 0', except for the points for which — 0' ~ 0, 7r, where the phase 5 
becomes relevant. A neutrino factory is expected to be sensitive to sin #13 > 1CT 4 [22] 
and to be able to see CP violation for phases of order 1 if sin #13 > .01 [23]. 

These plots can be understood analytically using the approximation for the CP 
asymmetry, eq.(30). When V L = 1, the CP asymmetry has a fairly simple expression 
in terms of low energy neutrino data and y 1 , the lightest eigenvalue of Y^Y^: 



%| Im ( m »i c i3 c i2 e i4> + m l 2 c i3 g i2 g W + ml 3 s 2 13 e 2lS 
8ttD I m Ul c 2 13 c 2 12 e-^ + m V2 cf 3 s? 2 + m U3 s 2 13 e~ 2iS ' 



where D = m 2 Ui c\ 3 c\ 2 + rn 2 2 c\ 3 s\ 2 + m l 3 sf 3 . For the mass hierarchy and the ranges 
of CHOOZ angles that we are using, it turns out that m U2 ^> m Ul , m V3 s 2 3 , so we can 
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Figure 1: "Phase overlaps" for the case Vl = 1, i.e. YJY,, diagonal. The upper left plot indicates 
the meaning of the distances to the different sides of the triangle. The rest of the triangles show 
density plots of the "phase overlaps" for random values of the phases and different CHOOZ angles: 
0.1 (upper right), 0.01 (lower left) and 0.001 (lower right). The darkest regions correspond to the 
largest density of points. 
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expand the denominator in e. Approximating #12, #23 ~ 7r /4, the result is: 



24 MM • (45) 




When the CHOOZ angle is much larger than ^m VY m 2 v Jm z Vi the first term in eq.(45) 
dominates, unless 25 — 0' is close to or tc. This condition is satisfied in particular 
when the lightest neutrino is very light, which is an interesting physical possibility. For 
the mass hierarchy that we have chosen, the condition above reads S13 ^> 0.01, which 
is satisfied when the CHOOZ angle is close the present experimental limit. Recall that 
#13 > .01 is required to detect 5 at a neutrino factory, so this limit would hold if CP 
violation is observed at the neutrino factory. The CP asymmetry in this case can be 
approximated by 

o..2 /„ \3 

sin(2<5-0'), (46) 

that does not depend on 0; only on 5 and <f)'. Furthermore, the dependence is such 
that one cannot conclude whether the " leptogenesis phase" is 5 or 0'. Instead, in 
this limit, the "leptogenesis phase" is the combination 25 — 0'. Comparing eq.(46) 
with the Fourier expansion, eq.(40), it follows that for most points, e ~ Csp. Hence 
e<£ ~ 0, e$ — e<f>' — Csft/2. (e$ ~ is a consequence of the fact that e depends on a 
combination of 5 and 0'.) Consequently, most points in the scatter plot, Fig. 1, upper 
right, are concentrated in the middle of the side corresponding to = 0. 

On the other hand, when the CHOOZ angle is very small, it is the second term in 
eq.(45) the one that dominates, as long as (0 — 0') is different from or n. In this 
limit, 

e~M I^! sin (0_0') . (47) 
4tt m V2 

The CP asymmetry only depends on and 0', and the "leptogenesis phase" is — 0'. 
As before, comparing eq.(47) with the expansion eq.(40), we conclude that for most 
points, e ~ C^y. Hence, e$ ~ 0, ~ ~ In consequence, most points in 

Fig.l, lower right, are concentrated in the middle of the side corresponding to 0$ = 0. 



Finally, for values of the CHOOZ angle ~ ym^m^/m^, both terms in eq.(45) have 
to be taken into account. In this case, we cannot say that there is a single "leptogenesis 
phase": both 5 — 20 and — 0' are "leptogenesis phases". Concerning the contributions 
from the phases 5, and 0' to the CP asymmetry, it is apparent from eq.(45) that in 
a generic point the three contributions are going to be comparable. To be precise: 

£5 = 

e 4> = 2 C W ( 48 ) 
e 4>' = 2^ s ^' ^H 1 ) 
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where, 



Css' — 



M 

3y 2 m 



m 



m 



'13 



sin(2<5 - 0') 



(/ 2 , 



"1 



47r m 



sin(0 — 0') 



;/ 2 



From these formulas, it is possible to understand the circular pattern that appears in 
Fig. 1, lower left, changing from triangular coordinates to Cartesian coordinates. We 
define the Cartesian axes setting the origin at the lower left vertex of the triangle, and 
we denote as x (y) the horizontal (vertical) axis. The change of variables read, 



2 

x = —= 

y = o] 



n 2 



(49) 



and using that ~ e$ + e^, we obtain, after some algebra, (x — -^) 2 + (y — |) 2 ~|, 
which is the equation of a circle centered in the barycentre of the triangle, with radius 
1/3. 



7 The general case 

In the general case the number of unknown parameters involved is rather large (five 
phases, two angles in the Vl matrix and the CHOOZ angle), so the analysis is much 
more intricate since many different limits arise. However, we will see that only a few 
limits are distinct and physically interesting; the rest correspond to small regions in the 
parameter space that could arise in particular models, but that we will not consider, 
following the same bottom-up spirit as in the rest of the paper. 

The different limits stem from the possible ways to expand the denominator in our 
approximate expression for the CP asymmetry 

3y 2 f E l W^ml ) , x 



The relevant elements in W for the calculation are 

"1 r r ...... r r fl 1 \ r /l 1 



±= 44 + e^ 2 44 (I + ^13) - ^ 13 4 (I - \e iS s 13 



-j= 44 - e^-44 (2 - 2 elSs ^) + (2 + 2 et5s 



13 



(51) 



1 .-,„„ r. r 1 



13 e c 12 c 13 i, 13 12 C 13 3 



where we have approximated cos #13 ~ 1 and we have assumed maximal solar and 
atmospheric mixings. It is apparent from these equations that different limits are going 
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to arise depending on the mixing angles in Vl and the CHOOZ angle. We find then 
an interesting interplay between leptogenesis and lepton flavour violation, induced by 
radiative corrections through Y^Y^ = V l \D\Vl- However, from the parametrization 
we have chosen for Vl, eq.(21), one realizes that the off-diagonal elements of Y^Y^ 
depend also on 9^, that does not play any role in the CP asymmetry generated in the 
decay of the lightest right-handed neutrino 4 . 

We obtain simple analytic expressions when the denominator of eq. (50) can 
be expanded in small parameters. In the whole of section 7, we assume \Wu\ > 
/m V2 , so we can neglect m Vl terms in the denominator. When IW13I 2 < 



\W'~ V " 



ni 



12 



ml /ml which will be the case if 6*f 2 , #13, 6*13 < -1, e can be approximated by 



3y? 



8?r W12 



-Im 



m. 



m 



V2 , 



m 



m 



V 2 , 



w? 2 



(52) 



On the other hand, when the mixing in Vl is large, in the sense that sin Q\ 2 or sin Q\ z 



is larger than ~ 0.1, then IW13I > \W- 



121 



*m 



1'2 



/m U3 , and the CP asymmetry reads 



3y! 



m,. 



8vr|W^i 3 | 2 \m 



Im 



v 3 . 



w? 



13 



(53) 



There is also an intermediate case, between these limits, where \Wi^\ 2 m u .Jm V2 < 
\Wi 2 \ 2 < \Wi 3 \ 2 ml J m 2 v . We do not discuss this, because m Vz jm V2 ~ 10 for the 



hierarchical LMA solution we consider. 
Let us analyze the two cases separately. 



7.1 \W n \ ~ IW12I » \W u \m v Jm V2 

The analysis for this case is parallel to the one we performed in the previous section for 
the case Vl = 1, where \Uu\ ~ \U\2\ \U\z\m U2i /m V2 . Using that m V2 ^> m Vl , m V3 SiSj, 
where is any of S13, sf 2 , s^ 3 , we can expand eq.(52), keeping the leading order terms 
in the expansion. The result is: 

\ 3 
Tft \ r 

[2s 2 13 sin(25 - 0') - (4) 2 sin(2y9 12 + 0') - (4) 2 sin(2^ 13 + 0') 

m V2 ) 

— 2\/2si3Si2 sm (^ - ^12 - 0') - 2"\/2si3Si3 sin(5 - </?i3 - 0') - 2sf 2 sf 3 sin(v?i2 + </?13 + 0') 

Obviously, in the limit sf 2 , sf 3 — > we recover eq.(45). In the case V L = 1 we found dif- 
ferent limits, depending on the value of the CHOOZ angle. Now, the role of the CHOOZ 

4 It is important, though, for the computation of the CP asymmetry generated in the decay of 
the heavier right-handed neutrinos, that could be relevant, or even dominant, in some scenarios 
(particularly in scenarios with non-thermal creation of right handed neutrinos). Research along this 
lines would be certainly interesting, since in this case lepton flavour violation could be intimately 
related with leptogenesis. 
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angle is played by the angles 9[ 2) #f 3 and the CHOOZ angle itself, and the results are 



"in- 



different depending on the values of these angles compared with ^m Vl m 2 U2 / 

• When any of the angles Si 3 , sf 2 or sf 3 is much larger than yj^n^m^/m^, the 
term proportional to (m U3 /m U2 ) 3 in eq.(54) dominates: 



2s 2 13 sin(25 - 0') - {s L l2 f sin(2<^ 12 + 0') - {s L l3 f sin(2^ 13 + 0') 



-2v / 2si 3 sf 2 sin(5 - </?i 2 - 0') - 2^513513 sin(5 - (p 13 - 0') - 2s[ 2 s[ 3 sin(</?i 2 + </?i 3 + 0') 

We recall here that this limit corresponds to the case where the lightest neutrino mass 
is very small. On the other hand, for the mass hierarchy that we are using as reference 
to present our numerical results, m Ul : m U2 : m U3 = 10~ 2 : 0.1 : 1, this limit arises when 
any of the angles is much larger than ~ 0.01, in particular, when the CHOOZ angle 
is close to the experimental bound and the relevant angles in Vj, are comparable to or 
smaller than the CHOOZ angle. 

When the three angles are comparable in size, we see that there are three "lepto- 
genesis phases": 25 — 0', 2tp 12 + 0' and 2<y? 13 + 0' (the arguments of the sines in the 
last three terms of eq.(55) are combinations of these). Notice that in this limit 0' is an 
important phase for leptogenesis, although it cannot be regarded as the "leptogenesis 
phase" , since the actual "leptogenesis phases" are combinations of 0' with other phases. 
However, an indication for a non- vanishing 0', coming for example from experiments 
on neutrinoless double beta decay, would provide an indication for leptogenesis. 

If there are two angles that are comparable, while the third is much smaller than 
the others, then there are two "leptogenesis phases" . To understand better the results 
for this limit, we analyze in some detail the case S13 ~ sf 2 ^> sf 3 . If s 23 is also small, 
this case would produce small rates for /i — > ey, as can be checked from eq.(21). The 
results for the other possibilities, S13 ~ sf 3 3> sf 2 and sf 2 ~ sf 3 3> S13, can be easily 
deduced from this analysis, making the appropriate substitutions. We have computed 
numerically the different contributions to the CP asymmetry for the choice of angles 
s 13 = sf 2 = 0.03, sf 3 = 0, the mass hierarchy m Ul : m„ 2 : m U3 = 10~ 2 : 0.1 : 1 and 
assigning random values, between and 2tt, to the phases. We obtain that for most 
of the parameter space, the only non-vanishing contributions to the CP asymmetry 
are e<j, e^ 12 and (</?i 3 does not play any role, because we have set s\ 3 to ). Since 
there are essentially only three contributions involved, a convenient way of presenting 
the results is using a triangular plot. In Fig.2, left, we explain how to interpret the 
distances to the different sides of the triangle, whereas in Fig.2, right, we show the 
numerical results of the calculation. We find that in general the three contributions 
are comparable, although the contribution from 0' is slightly larger than the other two. 
This can be understood from the analytical approximation, eq.(55), setting sf 3 = 0. 
The different contributions to the CP asymmetry are: 

e<5 = + gC^Vw 

e 4>' = 2^ 5< t>' + ^Vu) + gC^Vw (56) 
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Figure 2: The same as Fig.l, for the case where S13 = s{ 2 = 0.03, sf 3 = 0. The left plot indicates 
how to interpret the distances to the different sides of the triangle, and the right plot shows a density 
plot of the "phase overlaps" for random values of the phases. The darkest regions correspond to the 
largest density of points. 



where, 

Arc \m U2 ) 

C^ 12 - |f 3 (4) 2 sin(2^ 12 + 0') (57) 

that are in general comparable. 

Finally, if one of the angles dominates over the others, the conclusions are very 
similar as for the case Vl — 1, where the CP asymmetry received contributions from 
0' and 5. Here, the role of 5 is played by the phase corresponding to the angle that 
dominates (5 for s 13 , <p 12 for sf 2 , and (/? 13 for sf 3 ). In this case, e ~ C X( p', where x is the 
relevant angle among 5, <^ 12 and <y? 13 . On the other hand, the normalized contributions 
are O x ~ ~ 1/ \/2, while they are vanishing for the rest of the phases. For example, 
if sf 2 > si 3 , sfg, then e ~ and 0^ 12 ~ Op ~ 
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• When all the angles (S13, sf 2 and sf 3 ) are much smaller than yrn^rn^Jrn^, the 
term proportional to m Ul /m U2 dominates in eq.(54) and the CP asymmetry reads 

e „ M h<\ sin(0 _ ') . (58) 

47r ym^y 

In this limit, the results are identical as in the corresponding limit in the case Vl — 1, 
and there is a single "leptogenesis phase" ,0 — 0'. So, the normalized contributions to 
the CP asymmetry from and 0' are equal to 1/ v^2, while the contributions from the 
rest of the phases vanish. The numerical analysis yield a plot that is very similar to 
Fig.l, lower right, where the role of 0$ is played by either Vl2) O vl3 or Os- 

• Lastly, in the situations where the two terms in eq.(54) are comparable, the 
analysis is very involved, since in principle there are four independent "leptogenesis 
phases", namely, 25 — 0', 2</?i 2 — 0', 2</?i 3 — 0' and — 0'. So, the CP asymmetry 
receives contributions from the five phases, and in general they are comparable in size. 
Hence, it is very difficult to extract any general conclusion for this case. 



7.2 \W l2 \ 2 mJm Vs < \W U \ 2 

For simplicity, and since the number of phases and angles involved is rather large, 
we will set one of the angles in Vl equal to zero, say sf 3 = 0, so the phase </?i 3 
becomes irrelevant. Since sf 2 an d sf 3 appear in a similar way in the formulas, one can 
qualitatively derive the result when sf 3 is different from zero. With this choice, we 
are left with only two angles, the CHOOZ angle, S13, and one angle in Vj,, sf 2 . The 
limit we are studying in this section requires sf 2 larger than ~ 0.1. Then, using the 
experimental bound on the CHOOZ angle and that our phases are generically of order 
1, the denominator can be expanded as 



1 



13 



is L l2 ) 2 



(^1 +2^^513 COs(5 + y?i 2 )^ . 



(59) 



Hence, the CP asymmetry can be approximated by 



3yf (c 



2 („L \3 
12) 




^ -2 sin(2v?i2 + 0') + 2V2 (^j sin(v?i2 + 0') 




47T { S [ 2 f 

+S13 [2v^[sin(<J - y?i2 - 0') - 3 sin(5 + 3v?i2 + 0')] 

-4 [sin(5 - 0') - 3 sin(5 + 2cp 12 + 0')] + V2 (^j [2 sin(5 - ip 12 - 0') 



(60) 



'12 



+ sm(S + <p 12 - 0') - 3 sin(<5 + v?i2 + 0')] - 2 ( \ sin(5 - 0') 



This expression is rather cumbersome and it is difficult to extract information from 
it. It is not possible in general to identify the "leptogenesis phase", although it is 



21 



clear that leptogenesis depends mainly on 0' and </?i2, whereas the dependence on 5 is 
weaker. 

In Fig. 3 we show the numerical results for this case. As usual, we show a triangle 
with the meaning of the distances to the different sides (upper left plot), and density 
plots of the "phase overlaps" for the mass hierarchy m Ul : m U2 : m ua = 1CT 2 : 0.1 : 1, 
taking random values for the phases between and 2n. In the upper right (lower left) 
plot we show the results for tan#f 2 = 0.5 (1) and s i3 = 0.1. In both plots, the points 
are concentrated close to the base of the triangle (that corresponds to 0$ small), due 
to the small value of the CHOOZ angle. In the plot corresponding to tan#f 2 = 0.5 the 
points are concentrated around the center of the base, whereas for tan 6^ = 1, they are 
spread all over the base. This can be understood from the dependence of e on cot#f 2 . 
For tan#f 2 = 0.5, the terms with both </? 12 and 0' are the dominant ones, so C(f>'ipi 2 is 
the largest contribution to the CP asymmetry. On the other hand, when tan#^ 2 = 1 
these terms are comparable to the one with sin0', so e is dominated by C^^^ and C^. 
Depending on the value of 0' the points spread along the basis of the triangle. In the 
lower right plot we show the numerical results for tan#f 2 = 0.5 and s 13 = 0.01. The 
plot is similar to the one with S13 = 0.1 but with an even smaller value of Os- 



If neutrino masses are due to the seesaw mechanism, then the heavy right-handed 
neutrinos can generate a lepton asymmetry in the early Universe when they decay out 
of equilibrium, if they have CP violating couplings. Such complex couplings in the 
high energy parameters could induce three phases in the light neutrino sector, called 
0, 4>' and 5 (these are the phases that appear in the MNS matrix; see eqs. (7) and (8)). 
Upcoming experiments may be sensitive to two of these phases: the Dirac phase 5 could 
be measured at a neutrino factory, whereas the Majorana phase 0' might have some 
observable effects in neutrinoless double beta decay. In this paper, we are interested in 
the relative importance of the phases 0' and 5 for leptogenesis. 

To address this issue, we use a parametrization of the seesaw in terms of weak 
scale variables: the light neutrino masses, the MNS matrix, the eigenvalues of the 
neutrino Yukawa matrix, and a unitary matrix Vl- We assume a hierarchical light 
neutrino spectrum, with the lightest neutrino mass of order m^/100, and an MNS 
matrix that corresponds to the LAMSW solution to the solar neutrino problem. The 
matrix V L is related to the off-diagonal (lepton flavour violating) elements of the slepton 
mass matrix, and contains three phases, two of which (y?i2,y?i3) are relevant for our 
calculation. It is important to use a parametrization in terms of "physical" weak scale 
phases; this is discussed in section 5. 

In the parameter space we are interested in, we find a simple analytic approximation 
for the lepton asymmetry e±, produced in the decay of the lightest right-handed neutrino 



(see eq. (30)). In this equation, y x is the smallest eigenvalue of the neutrino Yukawa 



8 Summary and Discussion 




(61) 
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Figure 3: The same as fig.l for different situations where \Wi2\ 2 m U2 /m„ 3 < |Wi3| 2 . The upper 
left plot indicates the meaning of the distances to the different sides of the triangle. The rest of the 
triangles show density plots of the "phase overlaps" for random values of the phases, sf 3 = and 
tan = 0.5, Si 3 = 0.1 (upper right), tan#f 2 = 1, s 13 = 0.1 (lower left), and tan#f 2 = 0.5, s 13 = 0.01 
(lower right). The darkest regions correspond to the largest density of points. (See sect. 7.2 for 
details.) 
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matrix, and W is the unitary transformation from the basis where the vl mass matrix 
is diagonal to the basis where YjY^ is diagonal: W ln = [V L ]i m [U] mn where V L and U 
are defined in eqs. (7), (8) and (21). 

We are interested in the relative importance of the phases 0' and S for e. That is, 
we do not discuss whether we get e large enough, which is essentially controlled by real 
parameters, such as y\. We assume that the observed baryon asymmetry is produced 
in the out of equilibrium decay of z/r x and study how important could be 5 or 0' for 
the CP asymmetry. In other words, if we suppose that e is of the correct size, what 
fraction of e is due to S or 0'? 

We Fourier expand e on the five relevant phases of our low-energy parametrization 

e = A ikimn sin(j<5 + k(j> + 10 + rrup 12 + mp 13 ) , (62) 

j,k,l,m,n 

and then divide the sum into five components, one due to each phase. In e,y, which 
is the component due to S, we put all the terms from the Fourier expansion that are 
oc sin(j5). We define C$ a to be the sum of all the terms oc sin(j<5 + na), and divide it 
equally between e$ and e a — that is, we add \ ^2 a Cs a to e<5- We also add the terms 
oc sm(j5 + na+mfl), multiplied by 1/3, and so on (a and (3 are one of {0, 0', ipi 2 , fn})- 
This procedure is described in more detail in section 5, and the formula for e«5 can be 
found in eq.(37). Then we define a normalized "fraction of e due to 5" , which we call 
the overlap between the leptogenesis phase and 5, as 

O s = ^M= . (63) 

The magnitude of Os or Op depends on five phases and three unknown real parameters: 
the CHOOZ angle, 9 13 , and two angles from Vl, that is related to radiative decays. 
In the numerical calculation, we fix the hierarchy of light neutrino masses, the angles 
and assign random values (linearly distributed) to the phases between and 2n. The 
numerical results are shown in density plots in the Os — 0$ space. 

We present results for two representative cases. In section 6 we discuss the Vl = 1 
case, where the three relevant phases are 5, 0' and (the phases of the MNS matrix), 
and in Section 7 we allow for two non-zero angles in Vl- For the sake of clarity in 
the presentation, in Section 7 we analyze simplified scenarios to reduce the number of 
phases involved to three. Since the definition of overlap satisfies the identity J2 0\ = 1, 
a convenient way of showing the results is by using a triangular plot, where the distance 
to each side of the triangle corresponds to 0\. 

The Vl = 1 model should be a good approximation when the angles in Vl are 
smaller than the CHOOZ angle. It can be checked from eqs. (51) and (61) that the 
Dirac phases 5, tpyi and </? 13 appear in e multiplied by the sine of a real angle (#13, Q\ 2 
and respectively). If Q\ 2 and #f 3 are much smaller than the CHOOZ angle 13 , then 
</?i2 and </?i3 are less important for e than 5. 

In Section 7 we analyze the situation where there are angles in Vl larger than the 
CHOOZ angle. The angles of Vl are related to the branching ratios for £j — > ^7, 
as discussed after eq. (22) : BR(£j — > ^7) oc yfc|[Vj,]fc. ? -[Vi ( ]fcj| 2 . Current limits on 
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t — > /ry, r — > ej, and — > cy are satisfied if all angles 9^ < .1. However, present 
bounds and anticipated improvements on all three branching ratios can be satisfied if 
e.g. 0i 3 ,0% 3 — and Q\ 2 ~ 1. So it is phenomenologically possible to have at least 
one large angle. The associated phase could then be important for leptogenesis; this 
possibility is studied in that section. 

The approximate expression for e given in eq.(61), shows that for generic Wij, 
e oc Im(W^i 2 W / 1 * 3 ) 2 . That is, the terms proportional to m Ul can be neglected when 
Wu > Wn\Jm Vl /m V2 and W i3 > W\\^Jm Vx m 2 U2 / m\ z , as discussed in section 7. Recall 
that W 13 contains terms ~ sin6> 13 , sin#f 2 , sin#f 3 . We set sin#f 3 to zero, to ensure that 
^3 — * ^il constraints are satisfied, and because the functional dependence of e on sin Q\ 2 
and sin 9^ 3 is similar 5 . So the case studied in that section has three phases: 8, 0' and 

The CP asymmetry e usually depends on the interference between at least two 
phases. Our results can be divided into three representative cases, according to the 
neutrino masses and low energy mixing angles. The phases relevant to leptogenesis can 
be 1) the light majorana phases and 0', or 2) the neutrino factory phase S and 0', or 
3) 0' and phase(s) from the slepton mass matrix. The first two cases were discussed in 
section 6, and the third in section 7. We outline here the observational consequences 
of each case. 

The majorana phases and 0' are the relevant phases for leptogenesis (equivalently 
and Op are large) when 



«i3, »ij « 



\ 



ml 2 m vi 
ml 3 



, ¥ large) (64) 



where j = 1,2. For m ui ~ .01m u:i , 6 these conditions imply a CHOOZ angle #13 < .01, 
BR(t -> ei) < 1CT 10 and BR(n -> erf) < IQ~ a BR(t -> ^7). 7 So if this scenario is 
realised, then with forseeable sensitivities, r — > ^7 and /i — > could be observed, but 
9 13 and r — > e7 will not be seen. If was measured in 0i/2/3 decay, then sufficient 
conditions (eqn 64) could be determined: knowing m Ul fixes how small 6 13 and Q\- must 
be for their associated phases to be irrelevant. 

Notice that r — > 7x7 has little impact on the importance of <5 for leptogenesis. This 
is because it is mostly related to #23 > the angle from Vl which does not affect e in our 
p ar amet r isat ion . 

A second possibility, where 5 would be important for leptogenesis, arises if 



\ 



777,2 jji 

« i3 > iPs, 0# large) (65) 

m, 3 



5 this is because ^3 ~ [y T + v^)/^/2 

6 We use m Vl = .01m U3 to estimate numerical upper bounds on the angles. The parameter space 
where this scenario obtains shrinks down as m Vl decreases. 

7 In this discussion, we approximate (!) BR(£j -» t a ) ~ 10"' 1 yi\Vl kj V Lki \ 2 (^) (j^f) > so 
these estimates are for flavour only and should be taken with large crystals of salt. 
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This in the interesting case for the neutrino factory; it would arise for #13 > .01^/-^^— 

and 6fj < #13. For m Vl ~ .01m U3 , it corresponds to a CHOOZ angle accessible to the 
neutrino factory, r — > ey below anticipated sensitivities (BR(t — > e^y) < 1CT 9 ), and 
BR(fi — > e7) < 9\ z BR(t — > /ry). Case 2 (eq 65) corresponds to 0$ large, and therefore 
r — > e7 unobservably small. 

A final possibility is that the phases of the slepton mass matrix are more important 
for leptogenesis than 8. This corresponds to 



7 13, 



\ 



m 3 



«0£ , (0^,0^ large) (66) 



An example of this case would be 6>i 3 ~ .01 — so detectable at the neutrino factory, and 
sin^fg ~ — which would imply BR(t — > &y) ~ 10 -8 . This example is plotted 

in figure 3 (if {di 3 ,(pn} and {812,^12} are interchanged). If r — > cy is observed, then 
we are in case 3, Vlj > 0$ and the slepton phases are probably more important for 
leptogenesis than 8. Notice however, that lepton flavour violating processes could be 
small even in this case, where 8 is not important for leptogenesis. The 1-2 angle of 
Vl can be large, without inducing observable £j — > ^7 rates, because it induces lepton 
flavour violating slepton masses proportional to y\. Such simultaneously small lepton 
flavour violation and 0$ could be disfavoured by model-building, because it requires 

013 ^ 012- 

To summarise these three cases: we find that the neutrino factory phase 5 can 
be important for leptogenesis. When the CHOOZ angle is large enough to detect 
CP violation at neutrino factories (which requires #13 > .01), and the lepton flavour 
violating branching ratios are vanishingly small (#f 2 , #13 #13), then 8 contributes 
significantly to the leptogenesis phase ( | C*<5 1 2 > .3). This can be seen from our low- 
energy approximation to e, eq. (61): if Wi 3 is not too small, e depends on the phase 
difference between W 12 and W 13 , and if 9 13 > 9f 2 , #f 3 , the phase of Wi 3 is 5. Although 8 
appears always suppressed by the CHOOZ angle, which is small, it plays an important 
role, because W13 is multiplied by the largest neutrino mass m V3 . The phase 8 is 
unlikely to be important for leptogenesis if either the CHOOZ angle is small (#13 < .01 
— see fig 1), or if the r — > cy or fi — ■> cy branching ratios are large (6*f 2 , 9 12, > 6*13 — see 
fig 3). 

We find that the Majorana phase 0' is (almost) always important for leptogenesis. 
For instance, the fraction of e that is due to 0', 0<^, is significant in all the cases we 
have studied. Algebraically, the reason is that the main contribution to e in eq.(61) is 
generically proportional to Wi 2 , that is proportional to e 1 ^ , unless some cancellations 
occur. The contribution from any of the other phases can be suppressed by sending 
a small parameter to zero. The Majorana phase of m ui , <p, becomes unimportant as 
m Ul — > 0, and the three Dirac phases 8, ipu and ^13 multiply angles which are positively 
small (S13) or believed to be small (sf 2 > s i3)- The contribution of 0', on the other hand, 
is consistently significant. 

It is interesting to study how closely related are the Majorana phases of the light 
neutrinos to the Majorana phases of the heavy right-handed neutrinos. It is well 
known that when neutrinos have Majorana masses, there is CP violation even in the 
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two generation model. This suggests that the Majorana phases of the v R sector could 
be more important for leptogenesis than the Dirac phase, because they can contribute 
to the CP asymmetry e suppressed by mixing between only two of the u R , rather than 
mixing among the three z/r, as required for the Dirac phase. However, there is no 
symmetry-based distinction between Majorana phases and Dirac phases. The high 
scale Majorana phases are functions of all the weak scale parameters — the real ones 
as well as all the Dirac and Majorana phases. So the reason 0' is important is not 
that low energy Majorana phases determine the high scale Majorana phases. One can 
check from the formulae in section 4 that 0' is usually significant because it multiples 
a not-very-small mass, rather than a (possibly) small mixing angle. 

We do not find a simple correlation between the sign of low energy phases and the 
sign of the CP asymmetry e. Such a correlation would be interesting, and exists in 
certain models [45]. However, in our bottom- up approach, e is usually proportional 
to phase differences (e ~ sin(jo; + n/3)), as can be seen from eq. (61) and the various 
limiting cases discussed in sections 6 and 7. 

In summary, we have studied the relative importance of low energy phases for 
leptogenesis. Using a parametrization of the seesaw mechanism in terms of weak scale 
variables, we express the CP asymmetry produced in the decay of the lightest vr as a 
function of the "neutrino factory phase" S, the "neutrinoless double beta decay phase" 
0', and three other "physical" weak scale phases. We introduce a way of splitting e into 
contributions due to the different phases, 5, 0', etc. We assume that e is big enough to 
be responsible for the observed baryon asymmetry, and compare the relative size of the 
different contributions. We find that 0' is generically important for leptogenesis. The 
importance of 5 depends on the mixing angles of the slepton sector. If these are smaller 
than the CHOOZ angle, then 5 makes a significant contribution to leptogenesis. 
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Addendum 

After this work was completed, a related analysis appeared [51]. 

Appendix 

In the appendix we explain the numerical procedure that we have followed to compute 
the contributions to the CP asymmetry from the different phases. For the sake of clarity 
we will only present the procedure we followed for the case Vl — 1, where only the 
phases 5, and 0' were relevant. The extension to the general case is straight-forward. 
Our starting point to compute the contributions was the Fourier expansion of the 
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CP asymmetry 



e = A iki sm(j5 + k(f) + 1(f)') = C S + C <j) + C<p> + C &4> + C s< p> + + C SH> ' (67) 

j,k,l 

with C a , C a p and as in eqs.(34)-(36). From the periodicity of e it is apparent 

that 

Cs = jhfL d(t)d(t)lt 

C^j^J^dSdcf)' e (68) 
C,> = j±y 2 jy dtp e 



C54, = <¥ e-(C 5 + C#) 

C w = d<p e-(C s + C ) (69) 
1 

= j^—j J d5 e — (Ctj, + Op) 

Cs^' = e — (Cs + Cfj, + + Cstj, + Csft + C^i) . (70) 

These integrals can be computed numerically, thus giving the different contributions 
to the CP asymmetry. This avoids difficulties with the points where the approximation 
eq. (30) breaks down. However, it is also interesting to solve this integrals analytically, 
to cross-check the results we obtained in Section 6. The results for the double integrals 
are 

C 5 ~0 C ~O CV~0. (71) 

On the other hand, the results for the single integrals is more involved and depends on 
the particular point of the parameter space. These integrals can be computed using 
the residue theorem; the number of poles inside the unit circle depends on the values 
of the phases and other neutrino parameters, especially on the CHOOZ angle, hence 
the dependence of the result on the chosen parameters. However, some care must be 
exercised in using the residue theorem, because there can be poles in eq. (30) at points 
where the approximation breaks down. Such poles must be neglected. 

The results for the single integrals are different depending on the CHOOZ angle. 
We consider three possibilities: 

• When the CHOOZ angle is close to the experimental upper limit, or to be precise, 
when 

\ m vi c i3 c i2 ^ ~l~ m v A s i3 e 2l<5 | < m U2 c 2 3 s 2 2 

\m U2 c 2 13 s 2 12 e^' + m U3 s\ z e 2i5 \ > m Ul c 2 13 c 2 12 (72) 

vi ^13^12 c ^ " 6 ^2 u 13 d 12 c I ^ lih vz 1 
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the single integrals read 

^ " 87TD m \ m„ 2 c 2 3 s 2 2 eW + m U3 s 2 3 J " 
C^~0 , (73) 

where .D was defined after eq.(44). This result, coincide with eq.(46), that was obtained 
using a completely different method. 

• When the CHOOZ angle is very small, or when the conditions 

\ m ui C 13 C 12 e ^ "I" m vz s 13 e I *^ c 13 5 12 

l™m C ?3 S ?2 e* ' + s ?3 e2i<5 | < "V C ?3 C ?2 (74) 

l^fi c 13 c 12 "I" c 13 s 12 I ^ s 13 

are fulfilled, the results for the single integrals are 

CV ~ (75) 

n H t « c 2 3 c 2 2 e^ + m^ c? 3 sj 2 3y 2 . 

C^' ~ — — — — Im < — 1 — =— = ~l — n — 5 ~u } — -, sm(0 — <b) . 

w 8nD \m Ul c\ 3 c\ 2 e»* + m„ 2 c? 3 s 2 2 ' } 4vr m^ 2 vr r; 

This result is identical to the result obtained using series expansions in Section 6, 
eq.(47). 

• For intermediate values of the CHOOZ angle, it is usually the case that 

-.2 -o^^!) _i_ rT1 c 2 -j2i(5 1 /rr} ^2 c 2 

-'13 c 12 e t "^3 6 13 e I ^ Ul vi c 13 5 12 

c 2 3 s 2 2 e^' + m V3 s 2 3 e 2l<5 | > c 2 3 c\ 2 (76) 



TOj/j c 2 3 c 2 2 + 7?2j, 2 c 2 3 s 2 2 e l<?i | > m V3 s 2 3 
so the single integrals are 
Ou~0 



m»2 c i3 s i2 e^' + s 2 3 e 2i<5 

™3 2 2 i, 





8irD ' [ m Vl c 2 3 c 2 2 e i( t> + c 2 3 s 2 2 e^' J " 47r m 



3 

s 2 3 sin(25 - 0') 
-0') , (77) 



'^2 



that are identical to eq.(49). 
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